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Abstract 

This paper explores the use of laboratory closed-loop learning control to either fight or cooperate 
with decoherence in the optimal manipulation of quantum dynamics. Simulations of the processes 
are performed in a Lindblad formulation on multilevel quantum systems strongly interacting with 
the environment without spontaneous emission. When seeking a high control yield it is possible 
to find fields that successfully fight with decoherence while attaining a good quality yield. When 
seeking modest control yields, fields can be found which are optimally shaped to cooperate with 
decoherence and thereby drive the dynamics more efficiently. In the latter regime when the control 
field and the decoherence strength are both weak, a theoretical foundation is established to describe 
how they cooperate with each other. In general, the results indicate that the population transfer 
objectives can be effectively met by appropriately either fighting or cooperating with decoherence. 
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I. INTRODUCTION 



Control over quantum dynamics phenomena is the focus of many theoretical [1-3] and 
experimental [4, 5] studies. Various control scenarios exist, including optimal control[6, 7], 
coherent control[8], and STIRAP (Stimulated Raman Adiabatic Passage) control[9, 10]. In- 
creasing numbers of control experiments, including on complex systems[ll-18], employ 
closed-loop optimal control[19]. Many of the control systems explored theoretically are 
restricted to pure state dynamics or the dynamics of isolated quantum systems. In practice, 
control field noise [20-22] and decoherence[23, 24] inevitably will be present in the laboratory 
and the general expectation is that their involvement will be deleterious towards achieving 
control. Recent studies have investigated the effect of field noise and shown that con- 
trolled quantum dynamics can survive even intense field noise and possibly cooperate with 
it under special circumstances [25]. Decoherence of quantum dynamics in open systems, 
which often represents realistic situations, is a concern for control of atomic and molecular 
process, especially in condensed phases. Simulations have shown that it's possible to use 
closed-loop learning control to suppress the effect of quantum decoherence [26]. Some inves- 
tigations have been performed on laser control of population transfer in dissipative quantum 
system[27-30]. Control of decoherence and decay of quantum states in open systems has 
also been explored[31, 32]. Decoherence was shown to possibly be constructive in quantum 
dynamics [33, 34]. 

This paper considers the influence of decoherence (dissipation) upon the controlled dy- 
namics of population transfer with the decoherence induced by interaction with the en- 
vironment but spontaneous emission is not included. This regime arises, for example, in 
condensed phases where significant environmental interactions dominate the decoherence 
processes. The goal is to demonstrate that effective control of population transfer is possi- 
ble in the presence of decoherence. It is naturally found that decoherence is deleterious to 
achieving control if a high yield is desired, however we also find that good control solutions 
can be found that achieve satisfactory yields. If a low yield is acceptable, then it is shown 
that decoherence even can be helpful and the control field can cooperate with the decoher- 
ence. This paper will investigate these phenomena numerically and analytically to illustrate 
the issues, and this work compliments an analogous study considering the influence of field 
noise [25]. 
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The dynamical equations and control formulation is presented in Section II. Simulations 
of closed-loop management of dynamics with decoherence is given in Sec. III. Section IV 
develops an analytical formulation to describe how the control field and the decoherence 
cooperate with each other when they are both weak. Finally, a brief summary of the findings 
is presented in Section V. 



II. THE MODEL SYSTEM 



Realistic quantum systems in the laboratory often can not be fully described by a sim- 
ple, fully certain Hamiltonian. Instead, almost all real systems are affected by dissipation, 
decoherence or noise, which can not be totally eliminated. The recent successes of closed- 
loop optimization algorithms[19] operating in the laboratory demonstrate the capability of 
finding optimal, stable and robust solutions automatically [11-18], even for very complex 
systems. The present work aims to support the continuing experiments by investigating the 
principles and rules of control in the presence of decoherence. In keeping with this goal, a 
simple model system will be investigated. 

The effect of dissipation on controlled quantum dynamics will be explored in the context 
of population transfer in multilevel systems characterized by the dynamical equation of the 
reduced density matrix p s 



= -i[H - LiE(t),p s (t)} + 7^ { P s(t)} , (la) 
H = J2^\v)(v\, (lb) 

v 

where \v ) is an eigenstate of the controlled system Hamiltonian H , and e v is the associated 
-u-th field-free eigen-energy. The noise- free control field E(t) is taken to have the form 



E(t) =s(t)J2 A cos (uit + 6i) , (2a) 
i 

s it) = exp [- (f - T/2) 2 /2a 2 ] , (2b) 

where {ui} are the allowed resonant transition frequencies of the system, and off-resonant 
excitation are not included. The controls are the amplitudes {A{\ and phases {#;}, 
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and the control interacts with the system through the dipole operator p. In the laboratory- 
such control fields may be generated with programmable adaptive phase and amplitude 
femtosecond pulse shaping techniques [35, 36]. 

In Eq.(la) J 7 {ps{t)} is a functional which represents decoherence caused by the environ- 
ment and 7 is a positive coefficient which indicates the strength of the decoherence and will 
be varied in this paper to study the effect of dissipation on the control of quantum dynamics. 
Various equations have been developed[37-39] to describe the interaction between a system 
and the environment. In this paper we take the Lindblad form [3 7]: 

F{Ps(t)} = [ c {jn}Ps(t)C\ jn} - ^ Ps {t)C\ jn} C {jn) - \c\ jn} C {jn} ps{t)) , (3) 

{jn} V J 

Here C{j n } are bounded linear Lindblad operators acting on ps such that the norm of the 
operators, ^{j n } ^{jn}^0' n }> * s nnr t e - m t ne Lindblad equation, the Markovian and weak 
coupling approximations are made, and the normalization and positive definite nature 
of the reduced density matrix ps(t) is preserved. 

In the model system studied here, the Lindblad operators in Eq.(3) are expressed phe- 
nomeno logically as [40] 



C jn = ^T jn \j) (n\ , j, n = 0, 1, N - 1 , (4) 

where N is the number of levels and \j) are eigenstates of the isolated system Hamiltonian 
H . The coefficient Tj n , when multiplied by 7, specifies the transition rate between level j 
and n. These rates depend on how the system and environment interact and will be chosen 
arbitrarily for the purposes of the general analysis here. 

Inserting Eq.(4) into Eq.(3) produces the dynamical equation for the reduced density 
matrix of the multi-level system, 

dpivit) 



= -i{ei - e v )pu>(t) - i ^2 {VinPni>{t) ~ pin(t)Pni>) E(t) + ^J 7 {ps(t)} a , , (5a) 

n=0 

N-1 ^ N-1 



2 

n=0 n=0 



Closed-loop control simulations will be performed with model systems, but under the stan- 
dard cost function [20] realizable in laboratory circumstances: 
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J[E(t)] = \0[E(t), 1 ]-0 T \ 2 + aF, 



(6a) 
(6b) 



where Or is the target value (expressed as a percent yield) and 



0[E{t)n]=TrWf)O f ] 



(7) 



is the outcome produced by the trial field E(t) under decoherence strength 7, and F is the 
fluence of the control field whose contribution is weighted by the constant, a > 0. In the 
present work, Of — (tyf\ is a projection operator for the population in a target state 



III. NUMERICAL SIMULATIONS 

To demonstrate how quantum optimal control can either fight or cooperate with decoher- 
ence, we perform four simulations with simple model systems. The first three simulations 
use the single path system in Figure (la) which represents a model multi-level system[41] 
or a truncated nonlinear oscillator [42-44], while the last simulation uses the double-path 
system in Figure (lb). With the single path system, the first two models employ different 
decoherence coefficients associated with nearest neighbor transitions in Eq.(5b). The 
third model allows as well for next nearest neighbor transitions. The fourth model with 
the double-path system shows how cooperation between decoherence and the control field 
can influence the transition path taken by the dynamics. The decoherence term doesn't 
include spontaneous emission. All of the simulations employ closed-loop optimization with 
a Genetic Algorithm[45] (GA), in keeping with common laboratory practice. Equation (6a) 
is the cost function used to guide the GA determination of the control field. The model 
parameters below are chosen to be illustrative of the controlled physical phenomena, and 
similar behavior was found for other choices as well. 

A. Model 1 
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This model uses the simple 5-level system in Figure (la) with eigenstates i — 0, • • • ,4 
of the field free Hamiltonian H , having only nearest neighbor transitions with the frequen- 
cies c^oi = 1.511, uj\2 = 1.181, u 2 3 = 0.761, CU34 = 0.553 in rad-fs" 1 , and associated transition 
dipole moments /i i = 0.5855, fi 12 = 0.7079, /i 23 = 0.8352, /i 34 = 0.9281 in 10~ 30 C-m. The 
decoherence coefficients in Eq. (5) are also non-zero only between nearest neighbor levels: 

r i = r 10 = 0.0895, r 12 = r 21 = 0.1942, r 23 = r 32 = 0.1209, r 34 = r 43 = 0.2344. The 

target time is T = 200 fs, and the pulse width in Eq.(2b) is a = 30 fs. The control objective 
is to transfer population from the initially prepared ground state |0) to the final state |4), 
such that O = |4) (4|. In order to calibrate the strength of the decoherence, Table I shows 
the state population distributions when the system is only driven by various environmental 
interaction decoherence strengths. When 7 ^ 0.03 fs -1 , decoherence is very strong as it 
drives more than 30% of the population out of ground state at t — 200 fs. 

To examine the influence of decoherence when seeking an optimal control field, first 
consider a high target yield (i.e., perfect control with the target yield set at Ot = 100% 
in Eq.(6a)). The simulation results in Table II show that decoherence is always deleterious 
for achieving this target, but a good yield is still possible. In order to reveal the separate 
contributions of decoherence and the control field, Table II also shows the yield from the 
decoherence alone O [E{t) = 0,7] without the control field being present and the yield from 
the control field alone O [E(t),^ = 0] without decoherence being present. 

If we accept a low control yield outcome, very different control behavior is found in the 
presence of decoherence. The results from optimizing Eq.(6a) are shown in Table III with 
various levels of decoherence for a target yield of T = 5.0%. The results in Table III 
clearly indicate that there is cooperation between the control field and the decoherence. 
For example, when 7 = 0.03 fs -1 , the outcome of the control field and decoherence together 
(0[E(t), 7] = 4.94%) is much larger than outcome from control field alone (O [E(t), 7 = 0] = 
0.63%) or the decoherence alone (O [E(t) =0,7] = 0.67%), and even larger than their sum 
(1.30%). Table III also shows that the control process become more efficient with increasing 
decoherence strength (7) as the fluence of the control field is reduced while driving the 
dynamics to the same target. The control field spectra at different decoherence strengths are 
depicted in Fig. 2. The mechanism of cooperation between the control field and decoherence 
is not easy to identify from Fig.2. In the following simple model the decoherence coefficients 
Tij are specially selected to identify the mechanism of cooperation. 
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B. Model 2 



This model is similar with model 1. The only difference is that the transition rates 
in Eq.(5a) are chosen to increase monotonically: r i = Tio = 0.03495, r i2 = T 2 i = 0.1242, 
r 23 = T 32 = 0.3909, r 34 = r 43 = 0.6344. Table IV is the low target yield analog of Table III, 
but for model 2. Once again clear evidence is found for the optimal field cooperating with the 
decoherence. In order to deduce the mechanism of cooperation between the control field and 
the decoherence, Fig. 3 shows optimal control field spectra at different levels of decoherence. 
The control field consists of four peaks, corresponding to the four transitions between nearest 
neighbor levels. With increasing decoherence strength, 7, the peak intensities corresponding 
to transitions among the higher levels, become smaller, and successively disappear. This 
behavior shows that control field can be found to cooperate with decoherence and 
drive the dynamics more efficiently. As a result, the cooperation between the control 
field and environment becomes dramatic for 7 > 0.03 fs -1 . In this region the optimal 
field has essentially no amplitude at the 2-^3 and 3^4 transition frequencies such 
that the field alone produces a vanishingly small yield, yet with the environment present 
the cooperation is very efficient. If we chose the decoherence coefficients in Eq.(5a) 
decrease monotonically(not shown here), then with increasing decoherence strength, the 
peak intensities corresponding to transitions between the lower levels disappear successively. 
This behavior indicates a similar cooperation mechanism, but with a shift in the role of the 
various transitions. 



C. Model 3 



Model 3 is similar to model 1 except that two-quanta transitions are also allowed: oo 02 = 
2.692, c<Ji 3 = 1.942, cj 24 = 1.314 in rad-fs -1 , with transition dipole elements: /i 02 = —0.1079, 
A*i3 = —0.1823, /jl 2 4 = —0.2786 in 10~ 3 C-m. We also assume that there are environmentally 
induced two-quanta transitions: T 02 = 0.01099, T 13 = 0.1087, T 24 = 0.1346. Table V shows 
the outcome from model 3 with a low yield target Or = 10%. Once again, cooperation 
between decoherence and the control field is found, although to a lesser degree than in 
the earlier cases. The decrease in the degree of cooperation appears mainly to arise from 
the enhanced target value of Ot = 10%. As shown in model 1, eventually decoherence 
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has a deleterious effect when aiming towards a sufficiently high yield. However, under 
special circumstances, such as decoherence breaking the symmetry of the Hamiltonian H(t), 
decoherence could still have a beneficial role when seeking a high yield. 

D. Model 4 

Model 4 is the more complex two-path system in Figure 1(b). In this model, population 
can be transferred to the target state |4) along two separate pathways. The transition 
frequencies, decoherence coefficients and dipoles of the left path are the same as that of model 
1; the right path has the distinct transition frequencies: u>ov = 2.513, low = 1.346, = 
0.345,c^3/4 = 0.162 in rad-fs -1 , decoherence coefficients: T y = IVo = 0.1164, Fi> 2 ' = 
IVi' = 0.0885, IV3' = = 0.1557, IV4 = T43/ = 0.1280 and dipole elements: //or = 
0.6525, n vv = 0.7848, //^ = 0.9023, // 3 / 4 = 1.0322 in 10~ 30 C-m. For simplicity, only single 
quanta transitions are allowed along both paths. 

The results in Table VI show cooperation between the field and the decoherence in model 
4 for a low target yield of T = 5%. Figure 4 shows the power spectra of the three optimal 
control fields in Table VI. Panel (a) indicates that in the case of no decoherence the control 
field primarily drives the system along the right path. When decoherence is present in the 
left path (panel (b)), the optimal control field chooses to drive the system dynamics along 
the left path in order to cooperate with the decoherence. The lack of a peak at W34 in the 
latter case reflects the cooperation between the control field and the environment. Similarly, 
in panel (c), when decoherence is in the right path, the optimal field chooses to drive the 
system along the right path. These results clearly show that efficiently achieving the yield 
is the guiding principle dictating the nature of the control field and the mechanism in these 
cases. The cost function in Eq. (6) explicitly contains a fluence term which naturally guides 
the closed-loop control search towards efficient fields, and cooperating with decoherence is 
consistent with this goal. The common circumstance in the laboratory of having a fixed 
maximum fluence would produce similar behavior toward cooperating with decoherence 
when it's beneficial. 
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IV. FOUNDATIONS FOR COOPERATING WITH DECOHERENCE 



To illustrate the principle of how the control field can cooperate with decoherence consider 
excitation along a ladder (or chain) of nondegenerate transitions and energy levels, each 
linked only to its nearest neighbors as in models 1-2 in section III. One can also think of 
this system as a truncated nonlinear oscillator [42-44]. The N + 1 level system consists of 
an initially occupied ground state |0), N — 1 intermediate states \n), and a final target 
state | TV). The states are coupled with an external laser pulse having the nominal form in 
Eq.(2a). Both the control field and the decoherence are assumed to be weak in this section, 
corresponding to the low target yield cases in Section III. 

Just as the parameter 7 characterizes the strength of decoherence, it is convenient to 
introduce a parameter A to characterize the strength of control field in Eq.(la): 

dp sit) 



dt 



-i[H - p\E{t),p s {t)} + iF{p s {t)} . 



(8) 



In the analysis below we shall focus on the regime A — > and 7^0. Transforming the 
density matrix ps (t) to the interaction picture p (t) = e~ lHt>t ps (t) e lHot , or equivalently 



-i{ei-e v )t 



produces 



Piv = {ps) u < e 



^p(t) = [i\d (t) +1 F]p(t), 
C I (t)p(t)=[^(t)E(t),p(t)], 



(9) 

(10a) 
(10b) 



where pW (t) = e~ tHot pe tHot , or in matrix form py (t) = pkje~ lu)k ^ with uj k j = Sj — It 
is easy to show that the decoherence functional T defined in Eq.(5b) is not changed under 
the transformation of Eq.(9). The solution of Eq.(lOa) can be expressed formally as 



P ( T f) = exp+ 



(iA£j (t) + 7 J 7 ) dt 



P(0), 



(11) 



where exp + is the time-ordering exponential operator. In order to explore the scaling with 
7 and A we shall now introduce the first order Magnus approximation [46], which ignores the 
time ordering , such that 



P ( T f) ~ ex P 



(i\d (t) + 7 J 7 ) dt 



P(0) 



exp[(^ + 7 ^)T / ]p(0), 



(12a) 
(12b) 
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where 8 is a time-independent functional acting as 

£ P =[W,p] (13) 

with the matrix W defined as 

W=^ (t) E (t) dt, (14a) 

1 f Jo 

W kj = i- / ' m e- l ^E (t) dt » n kj ^^-. (14b) 
Here e (a;) denotes the spectrum of the control pulse: 

/oo 
E(t)e-^dt = J> w,) e*« + 5 (w + w,) e"*'] , (15) 

with g (uj) being the Fourier transform of the shape function s(t) in Eq.(2b). The first order 
Magnus approximation suffices in the "impulse limit" of a sufficiently short pulse[47]. 

It's convenient to introduce double bracket notation[48] where the ket \jk)) denotes the 
Liouville space vector representing the Hilbert space operator \j) (k\ and the bra ((jk\ as 
the Hermitian conjug ate to \jk)): ((jk \ = (\jk))^ = \k) From Eqs.(7) and (12b), the 
outcome from applying the control field can be expanded and rewritten as 

0[E{t),i\ = ((NN\exp[(i\S + 1 F)T f ] |00)) (16a) 

oo 

= (( NN i ( iX£ + pk ( T f) i °°» > ( 16b ) 

where the functional P is defined by: 

P [/](*)= / f(ti)dti (17a) 
Jo 

P" [/] (*) = I dh / 1 dt 2 --- [ dt k f (t k ) . (17b) 
Jo Jo Jo 

j-k 

It is easy to verify that P [/] (t) = y for the identity function f(t) = 1. 



A. Dynamics Driven by the Control Field Alone 

When there is no decoherence, the yield from the control field is 

OO rpk 

O[E(t),0]=J2 1 r((NN\( i £) k \00))\ k . (18) 
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The system yield can also be described in Hilbert space by 
O[E(t),0] = |(iV|e lH/T ' A |0>| 2 



OO rpk 



fc=0 



(19) 



For a ladder system, W is a tridiagonal matrix. When A — > 0, we only need to keep the first 
non-zero terms of Eq.(19), 

2 



O[E(t),0] 



rpN 

-j^(N\ (iW) N \0) 



rp2N 

\2N _ f rp2 \2N 



M 2 0N ' 



(20) 



where the coherent transition elements T mn are defined as 



T mn = \(m\W n - m \n)\ = J] |%_!| = [[^ 

fc=m+l 
k (^fc,fc-l)| 



k=n 



&k — \Wk,k-l\ — A*fe,Jb— l- 



/ 



(21a) 
(21b) 



By comparing Eqs.(18) and (19), we know that the first non-zero term in Eq.(18) is k = 2N 



2N 



T 

O[E(t),0]~^((NN\(i£) 2N \00))\™, 



(22) 



and the following equality holds: 



((NN\£™\00)) = C? N T* N . 



N rp2 



The extension of the above equality to additional matrix elements is straightforward, 



(23) 



{{n + m? n + m | {l£ fm | nn)) = CZT 2 



n+m' 



(24) 



Finally, the lowest order non-zero outcome from applying the control field without decoher- 
ence is 



O [E (t) , 0] = A 



-^2N 
2N ± f rp2 

N\ 2 0N ' 



(25) 



B. Dynamics Driven by Environmental Decoherence Alone 



For simplicity, the decoherence is assumed to only induce transitions between nearest 
neighbor levels in the system, so that the decoherence matrix T is tridiagonal with elements 



(26) 
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and T plays a similar role as the dipole matrix. Decoherence can also drive the system from 
the initial state to the target state, at least to some degree, when there is no external field. 
The outcome due to environmental decoherence without the control field is 

OO rpk 

O[0, 7 ] = ^-^(<A^ fe |00>>7 fe - (27) 
fe=i 

It is easy to obtain the following decoherence transition elements, 

TO 

((n + m,n + m\F m \nn)} = Y[-f n+k , (28) 

k=i 

where 7& is the fc-th decoherence-induced transition rate between the nearest neighbor levels, 

ik = r fc+1)fc . (29) 

If the decoherence coupling is weak, 7^0, only the lowest non-zero term of Eq.(27) is 
needed, 

rpN 

O [E (t) = 0, 7] « ^ ((NN I T N 1 00)) 7^ (30a) 

rpN N 

' k=l 

which corresponds to a transition path 1 00)) — > 1 11)) — > • • • — > \NN)). 



C. Cooperation between the Control Field and Decoherence 

Consider the case where the control field and environmental decoherence are present 
simultaneously, so that the outcome is described by Eq.(12b). If the control field and 
decoherence are both weak, the perturbation approximation can be used again. Strong 
cooperation between the control field (Eq.(25)) and decoherence (Eq.(30b)) is expected 
when their contributions to the outcome are the same order, which corresponds to 

7 = A 2 . (31) 

Then the perturbative solution p (t) has the property 

Pmn (t) ~ \ m+n . (32) 
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p(Tf) in Eq.(12b) is the solution of following equation 

dp(t) 



. }i ~(i\£ + \ 2 F)p(t), (33) 

expressed explicitly as 

w i\ {W h i-ipi-^ v - pi,v-iW V -i, v ) + \ 2 8u / rij- 1 p l -. 1 j- 1 , (34) 

where the terms higher than \ l+v are neglected. The last term of Eq.(34) represents the 
effect of decoherence which only drives transitions of the type |Z — 1,1 — 1)) — > \ll)). If a 
high control yield is expected and the perturbation theory approximation is invalid, all of 
the decoherence terms in Eq.(5b) have to be considered, where some terms will interact 
destructively with the control field induced dynamics, 

{Fp)ii> = ^w^i,i-iPi-i,i-i + ^ii'^i,i+iPi+i,i+i ~ 2 + ^1+1,1 + + r«'+i,z') Piv- (35) 

Here it's the first term that cooperates with the control field. 
The outcome from both the control field and decoherence is 

O [E (t) , 7] = ((NN I exp [(MS + \ 2 F) T f ] 1 00)) (36a) 

00 

= E « NN \ i XtSP ( T /) + a2 ^ p ( T /)] fc • ( 36b ) 



k=i 

v2N 



Expanding the above equation and keeping only terms of order A , we get the perturbation 
approximation for the outcome, 

N 

O \ E (*) , 7] ~ A 2Af E E A( ni ,n 2 ,--- ,n m ), (37) 

m=0 0<n 1 <- -<n m <N 

with A(n 1 ,n 2 , • ■ • , n m ) given by 
A (m, n 2 , ■ ■ ■ , n rn ) = P 2N ~ m (7» ((NN | (iSf N ~ n ^ \ n m n rn )) 

m 

x Yl (( n kn k I (iS) 2 ^-'^- 1 -^ I + 1, n fc _i + 1)) ((n fc _in fe _i | J" | + 1, ra fc _i + 1)) , 



fc=i 



(38) 



describing the contribution of the paths in which the decoherence operator T drives the 
system from state \rik — 1,^ — 1)) to |nfc,n fc )), and the control field operator £ drives the 
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system from state | n k + 1, n k + 1)) to \n k +in k+ i)). Substituting Eqs.(24) and (28) into 
Eq.(38), we have 

rp2N—m m 

A K n m ) = j^^C"-^ (39a) 

rp2N—m m 

_ rp2 f f-iN-Tim-1 TT ^ n fc /-)"fc-"fc-l-l l"?OM 

- J 07V / 2Ar _x,^2(7V-n m -l) 11^2 °2(n fc -n fc _i-l)- l dyD J 
V /• fc=1 n fc 

If each term in Eq.(2a) is only resonant with a single transition, 

u n « e n+ i - e n , (40) 

using the RWA[41] (Rotating Wave Approximation) we can drop all non- resonant terms in 
Eq.(lOa), 

fi {I) (t)E(t) ^ s(t)H F . (41) 
The time-independent Hp matrix elements are 

( H F)kj = VjAj5 kJ+1 + fi k A k 5 k+hj . (42) 

Here, /i k = /ik,k-i denotes the k-th dipole element between the nearest levels. The operator 
Ci {t) commutes with itself at different times implying that the Magnus approximation is 
valid in the present circumstance under the RWA, and it's easy to get Eq.(39b) with Q k 
being 

tt k = H k A k -^, (43) 
1 f 

defining the effective pulse duration as 

f T 



I ' s{t) dt. (44) 
Jo 



It is also easy to show that the outcome of Eq. (37) can be decomposed into a linear combina- 
tion of contributions from the field component intensity A? and the decoherence transition 
coefficients jj, 

0[E(t) )1 } = A)F l]+1] F 2v (45) 

where F X j and F 2 j are functions independent of Aj and jj, respectively. There is clearly 
cooperation between Aj (a coherently driven transition) and 7j (a decoherently driven tran- 
sition). For example, the outcome for a two level system is 

Pll (T / ) = (/i 1 A 1 ) 2 T e 2 + 7iT / , (46) 
14 



and the outcome for a three level system is 

P22 (T f ) = \ {^A^ 2 A 2 f T' e + i faAi) 2 l2T 2 e T f + l - (pt 2 A 2 ) 2 7l T e 2 7> + ^T} 2 . (47) 

The objective cost function in Eq.(6a) can be written in terms of the contributions from 
each specific control field intensity A 2 , 

J (A,) = (A)F l3 + l3 F 23 -0 T f + a A) + a £ A\ . (48) 

Mi 

Here, and F 2 j are independent of both Aj and 7^. It is easy to show that the condition 
for optimizing Eq.(48) with respect to Aj is 

4^' + ^ = T -^-, (49) 

The latter equation indicates that the contribution from decoherence can beneficially act to 
decrease the required amplitude of the optimal control field to attain the same yield. 

V. CONCLUSIONS 

Various impacts of decoherence upon quantum control are explored in this paper. Nu- 
merical simulations from several cases indicate that control fields can be found that either 
cooperate with or fight against decoherence, depending on the circumstances. Two extreme 
cases of high and low target yield are chosen to illustrate distinct control behavior in the 
presence of decoherence: (a) the control field can fight against decoherence effectively when 
a high yield is desired, while (b) in the case of a low target yield, the control field can even 
cooperate with decoherence to drive the dynamics while minimizing the control fluence. 

Four models were studied and a first-order perturbation theory of the weak control 
and decoherence limitis presented in this paper. Models 1 and 2 considered the control of 
the same system in different decoherence environments. Both cases show clear cooperation 
between the control field and the decoherence driven dynamics for low target yield. The 
detailed mechanistic role of decoherence can be subtle as indicated in model 1, although 
the final physical impact may be simple to understand. Model 2 clearly identified the role 
of decoherence by setting up a special structured interaction between the system and the 
environment. Although the conclusions in this paper are based on simple physical models 
with the Lindblad equation, similar behavior is expected for other systems and formulations 
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of decoherence. The findings in this work are parallel with an analogous study[25] on the 
role of control noise. These work aim to provide a better foundation to understand the 
physical processes at work when using closed-loop optimization in the laboratory, even in 
cases of significant noise and decoherence. 
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Table I. Population distribution of the single path model 1 when the system is only driven 
by interaction with the environment. 









Population in the state (%) 




7 "(fs- 1 ) 





1 


2 


3 


4 


0.05 


52.8 


25.5 


14.6 


4.64 


z.oy 


0.03 


65.0 


22.7 


9.65 


1.98 


0.67 


0.01 


84.8 


12.7 


2.25 


0.17 


0.02 


0.00 


100 














a Strength of decoherence, refer to Eqs.(l) and Eq.(5) 






Table II. Optimal control fields fighting against decoherence with model 1 for the highest 


possible objective yield of Or = 


100%. 








7 (fs" 1 ) O [E(t) 


= 0, 7 f(%) 0[E"(t),T = 


= 0] fe (%) O [E°*>(t), 7 ] c (%) O [E° p (t), 7 ] e (%) F d 


0.05 


2.39 


97.58 


34.90 




34.58 0.071 


0.03 


0.67 


98.43 


46.61 




46.42 0.067 


0.01 


0.02 


98.65 


72.90 




72.79 0.066 


0.00 


0.00 


98.53 


98.53 




98.53 0.064 



a Yield from decoherence alone without a control field 

b Yield arising from the control field without decoherence, but the control field is deter- 



mined in the presence of decoherence at the specified value of 7. 

c Yield from the the optimal control field in the presence of decoherence, Eq.(7). 
d Fluence of the control field. 

e Yield from decoherence and the field optimized with zero decoherence. 
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Table III. Yields from optimal control fields with different levels of decoherence in the 
case of model 1 for a low objective yield of Q T = 5.0%. 



7 (fs" 1 ) 


O[E(t) = 0, 1 ] a (%) 


O[£°P(t), 7 = 0] b (%) 


0[£^),7] C (%) 


p d 


0.05 


2.39 


2xl0" 4 


4.92 


4.08xl0~ 3 


0.03 


0.67 


0.63 


4.94 


8.17xl0~ 3 


0.01 


0.02 


3.01 


4.96 


1.23xl0- 2 


0.00 


0.0 


4.96 


4.96 


1.39xl0~ 2 



a,b,c,d refer tQ Table JT. 



Table IV. Optimal control with decoherence for model 2 with a low objective yield of 
Q T = 5.0%. 



7 (fs- 1 ) 


o[£(t) = o, 7 n%) 


O [£<*(*), 7 = 


0] 6 (%) 


0[^(t), 7 ] c (%) 


P d 




0.05 


2.2 


7.69x10- 


•12 


4.98 


2.07x10" 


-3 


0.03 


0.71 


2.51x10" 


-9 


5.08 


5.94x10" 


-3 


0.01 


0.03 


0.52 




4.90 


1.23x10" 


-2 


0.00 


0.0 


4.96 




4.96 


1.40x10" 


-2 



a,b,c,d refer tQ Table JT. 
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Table V. Optimal control with decoherence for model 3 with a low objective yield of 
Q T = 10,0%, 

7 (fs" 1 ) 0[E(t) = 0,7] a (%) O[£°P(t), 7 = 0f(%) 0[£°f(t), 7 ] c (%) F d 



0.05 5.24 3.49 9.92 1.03xl0~ 2 

0.03 2.08 4.73 10.00 1.17xl0- 2 

0.01 0.19 10.06 10.00 1.84xl0~ 2 

0.00 0.0 10.00 10.00 1.70xl0- 2 



a,b,c,d refer tQ Table JJ 

Table VI. Yield attained from the optimal field for model 4 with the low objective yield 
of Q T = 5%. 

7= (7L,7«) e (fs^) (0[E(t) =0, 1 ]) N "(%) O{E°?(t), 1 = 0} b (%) (O [E°*(t), 7 ]) c (%) F d 



0.00,0,00 


0.00 


4.99 


4.99 


1.18xl0~ 2 


0.04,0.00 


1.42 


1.34xl0" 4 


4.92 


6.20xl0" 3 


0.00,0.04 


0.85 


0.33 


4.95 


5.98xl0" 3 



a,b,c,d refer tQ Table JJ 

6 ILiIr'- Decoherence strength in the left and right paths, respectively; 
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FIG. 1: Two multilevel quantum systems used to investigate the impact of decoherence on optimally 
controlled dynamics, (a) The five level single path system used for simulations with models 1-3 
in Sections III.A-III.C, respectively. Models 1 and 2 only allow for single quanta transitions (solid 
arrows) while model 3 allows both single and double quanta transitions (dashed arrows), (b) The 
double-path system for model 4 in Section 3.D. 
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FIG. 2: Power spectra of the control fields for model 1 aiming at a low yield of Ot = 5.0% (Table 
III). 7 indicates the strength of decoherence (Eqs.(l) and (5)). The control field intensity generally 
decreases with the increasing decoherence strength reflecting cooperative effects. 
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FIG. 3: The same situation as Figure 2, but for model 2. Here the peaks corresponding to 
transitions involving higher quantum states disappear successively in keeping with the ability of 
decoherence to aid in meeting the desired physical goal. 
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FIG. 4: Power spectra of the control fields for model 4 (see panel (b) of Figure 1). The frequencies 
along the right path in Figure 1 are denoted with a prime to distinguish them from the left path, 
(a) No decoherence, (b) decoherence in the left path, (c) decoherence in the right path. In these 
low target yield cases the optimal field cooperates with the decoherence directing the dynamics to 
follow the associated path. 
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